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Abstract. Graphene has opened new avenues of research in quantum transport, with 
potential applications for coherent electronics. Coherent transport depends sensitively 
on scattering from microscopic disorder present in graphene samples: electron waves 
travehng along different paths interfere, changing the total conductance. Weak 
localization is produced by the coherent backscattering of waves, while universal 
conductance fluctuations are created by summing over all paths. In this work, we 
obtain conductance images of weak localization with a liquid-He-cooled scanning probe 
microscope, by using the tip to create a movable scatterer in a graphene device. This 
technique allows us to investigate coherent transport with a probe of size comparable 
to the electron wavelength. Images of magnetoconductance vs. tip position map the 
effects of disorder by moving a single scatterer, revealing how electron interference 
is modified by the tip perturbation. The weak localization dip in conductivity at 
B = is obtained by averaging magnetoconductance traces at different positions of 
the tip-created scatterer. The width ABwl of the dip yields an estimate of the electron 
coherence length at fixed charge density. This "scanning scatterer" method provides 
a new way of investigating coherent transport in graphene by directly perturbing the 
disorder configuration that creates these interferometric effects. 
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1. Introduction 

Graphene, a single atomic layer of carbon atoms in a hexagonal lattice (figure [TK), shows 
quantum phenomena pL] from the coherent fiow of electron waves. These include the 
quantum Hall effect ^ , the Josephson eflFect [3] and weak localization [4] . Graphene has 
an unusual band structure created by the interaction of electrons on the two sublattices, 
labeled A and B in (figure [IK): the conduction and valence bands are conical, and they 
meet at a point in fc-space, the Dirac point, like the band structure for a massless 
relativist ic particle. The ability to tune the charge density in graphene from positive 
values for electrons through zero to negative densities for holes, with no energy gap, 
means that the Fermi wavelength can be comparable to the mean free path, and the 
charge density can be broken up into puddles of electrons and holes by disorder. 

Coherent electron waves traveling through disordered conductors show weak 
localization j5|-0| and universal conduction fiuctuations (UCF) fs-lO]; both effects 



are created by the interference of electron waves. Weak localization results from the 
coherent backscattering of waves from a disordered potential (figure [lb) , while UCF are 
created by the interference of waves traveling along all possible paths (figure [l]!) . The 
set of paths that interfere is bounded by the difltusive coherent length ^ L^^. Because 
of the unique band structure of graphene, unique behavior for weak localization and 



UCF has be en p redicted jlll-115[|. and experiments are beginning to shed light on these 
issues BIliQ- 

In this work, we use a liquid-He-cooled scanning probe microscope to probe weak 
localization in graphene by mapping the effect of a single SPM-tip-created scatterer 
on coherent electron transport. While much can be learned from bulk transport 
measurements, a nanoscale probe that perturbs the system on the same size scale as the 
disorder potential and the electron wavelength provides important new information ^ 



30|]. By scanning the location of the scatterer created by the SPM tip, we spatially map 
the change AG in conductance G caused by the tip, producing a conductance image 
G(n, B) that represents a "fingerprint" of the intrinsic disorder, at a particular density 
n and magnetic field B. The magnitude of the conductance fiuctuations is ^ e^//i, their 
lateral size is ^ 10s of nm, comparable to the Fermi wavelength, and the images repeat, 
as found for UCF conductance images of graphene in zero applied magnetic field [16]. 
The conductance images change with magnetic field, and become uncorrelated when 
AB is larger than a characteristic correlation field ABc. Weak localization is observed 
as a dip in the magnetoconductance G at 5 = 0, with a characteristic width ABwl- 
We study weak localization by averaging over either a range of back gate voltages Vg to 
vary the charge density, or over a set of different tip positions at a fixed Vg. We find 
predicted agreement between the values of ABwl for weak localization and ABc for 
UCF and, as discussed below. Both ABc and ABwl increase as the charge density n is 
reduced, and display a maximum at the Dirac point. 
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Figure 1. (a) The hexagonal graphene lattice with carbon atoms in the two sublattices 
labeled A (red) and B (blue), (b) Schematic diagram of graphene's band structure 
showing how the two pairs of conical conduction and valence bands each meet at a 
point in /c-space - the Dirac point. The states are filled up to the Fermi energy Ejr. (c) 
Illustration of a pair of time-reversed backscattered trajectories that interfere to cause 
weak localization. The diffusive coherence length bounds the set of trajectories 
that can interfere, (d) A pair of forward scattering trajectories whose interference 
contributes to universal conductance fluctuations (UCF). 



2. Theory 

Graphene has an unusual bandstructure. A graphene sheet is composed of a single 
layer of carbon atoms in a hexagonal lattice, with two atoms in each unit cell forming 
sublattices labeled A and B in figure [IK- The band structure has two valleys: each 
valley has a conduction and a valence band that meet at a single point in the Brillouin 
zone, known as the Dirac point, with no energy gap (figure [lb) . The electron and hole 
energies increase linearly with momentum |k| near their meeting point, and they are 
nominally isotropic, yielding pairs of conical bands that are symmetric about the Dirac 
point. Because of the two atoms in each unit cell, there is an additional degeneracy 
in the electronic states of graphene known as chirality, with the eigenstates of opposite 
chirality arising from the two equivalent sublattices. Scattering between the two valleys 
in graphene and the effects of chirality must be considered in understanding weak 
localization in graphene, as discussed below. 

The phenomenon of weak localization can be understood by considering interference 
between time-reversed paths of backscattered electrons in a disordered material jsl, IgI . 
Figure [lb illustrates a pair of time-reversed paths that contribute to backscattering: both 
paths have an incoming wavevector kin, and an outgoing wavevector kout = —kin- The 
only difference between them is that one traverses the loop in a clockwise direction, and 
the other traverses the loop in a counterclockwise direction. As an electron travels along 
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each path, it accumulates phase according to the Dirac equation in the case of graphene, 
or the Schrodinger equation in more conventional materials. This phase is controlled by 
the integral of the electron momentum and the magnetic potential along the path, as well 
as the geometric phase acquired in graphene. In conventional materials such as metals 
and semiconductors with no spin-orbit coupling, the geometric phase difference between 
these paths around a closed loop is zero, yielding constructive interference at 5 = 
that enhances backscattering, causing a net reduction of the conductance (hence, weak 
localization). In graphene, the accumulation of geometric phase depends on whether 
elastic scattering breaks the chiral symmetry of the electron state or scatters electrons 
between valleys; these processes lead to both constructive and destructive contributions 
to the interference, as discussed below. 

An applied magnetic field acts to destroy weak localization by altering the phase of 
electrons traveling in opposite directions along time-reversed paths. In a perpendicular 
magnetic field S, a phase shift Scj) = ±27tAB /^q is picked up in both paths in figure [lb 
given by the magnetic fiux penetrating the effective area A enclosed by the diffusive 
loop, where = h/e is the magnetic fiux quantum - the sign depends on the direction 
of propagation around the loop. As the magnetic field is increased, the interference 
present at 5 = is destroyed, allowing one to observe weak localization experimentally 
as a dip in magnetoconductance at 5 = 0, with width ABwl determined by the average 
effective area A enclosed by pairs of backscattered paths. 

In the discussion above, we have assumed that the transport is entirely coherent, 
whereas realistically, inelastic scattering also occurs that randomizes the phase of 
an electron wave. If inelastic scattering occurs at a rate l/r^^ and we assume an 
elastic scattering rate I/t^ ^ l/'^</>5 then we can define a diffusive coherence length 
L(f) = {DTcf^y^'^^ where D is the electron diffusion constant; L^^ is the average distance 
over which an electron remains coherent. The coherence length L^^ sets an upper limit 
to the size of diffusive loops that give rise to weak localization (see figures [lb and [Hi), 
resulting in a characteristic magnetic field 

AB„^. « ^ (1) 

required to destroy weak localization, by producing a 27r phase shift between time- 
reversed paths in an area L^. 

The conductivity correction AgwL{B) due to weak localization is calculated by 
using the diffusion equation for coherent electron transport. In a two-dimensional 
conventional metal or semiconductor the change Ag in conductivity g vs. magnetic 
field S is 

where 7 is the total spin and valley degeneracy, and F{z) = ln{z) + '0(1/2 + 1/^), 
with '0(x) the digamma function. Equation [2] is calculated by considering the diffusion 
equation in the limit where the sample dimensions L^W ^ L^^, and the elastic mean 



Probing weak localization in graphene with a movable scatterer 



5 



free path 1^ <C Is- In graphene, equation [2] becomes 



12 



4eB 




(3) 

where = DL~^ is the rate of chirahty-breaking scattering within a vaUey, and 
T^^ = DL^^ is the rate of scattering between valleys. 

Because equation [3] has both positive and negative terms, the conductivity 
correction in graphene may be positive or negative at 5 = depending on the relative 
rates of the different scattering processes. A positive, weak anti-localization correction 
to the conductivity has been observed in epitaxially-grown graphene samples |20|. In 
exfoliated graphene samples, however, such as the ones considered here, the negative 
weak localization correction is typically dominant at low temperatures jsl. H 21 



implying the existence of a significant source of chirality-breaking scattering, possibly 
from lattice defects in the graphene. We observe a negative weak localization 
conductivity correction in all samples studied, and obtain qualitatively good fits to 
equation [3] with fit parameters similar to those in Ref. However, equation [3] is 
derived for L^W ^ L^^ which is not the case for our samples. In samples A and 
B (see figure [3b), W ^ L^fy^ and in sample C, we have W ^ L^fy near the leads and 
W ^ Lcfy in the central region. For simplicity, we will use equation [1] to estimate 
from measurements of ABwl below. 

In addition to weak localization, the interference of electron waves traveling 
along diflFerent paths creates universal conductance fiuctuations (UCF) with standard 
deviation SG ^ /h independent of sample size or the degree of disorder This 
coherent correction to the conductivity depends sensitively on the positions of the 
scatterers. We have recently shown jld] that moving a single scatterer created by an 
SPM tip by a distance comparable to the electron wavelength is sufficient to cause 
the full range of conductance ffuct nations in zero magnetic field, in agreement with 
theoretical predictions 



31, 32 



For UCF at finite B, the change A5c in magnetic field needed to reduce the 
correlation of the conductance by a factor of two is |[33|] 

^ A. (4) 

Note that ASc = 2 ABwl because weak localization arises from the phase difference 
accumulated along two counter-propagating paths around a diffusive loop, whereas 
UCF comes from a single loop. We compare these theoretical predictions with our 
experimental results below. 
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Disordered graphene 



Figure 2. Schematic diagram (to scale) of the charged scanning probe microscope tip 
creating an additional scatterer (red) in graphene. The shaded blue color scale in the 
graphene layer represents the disorder potential. 



3. Experimental Methods 

In this work, we use a liquid-He-cooled SPM tip to create a movable scatterer in a 
mesoscopic graphene Hall bar. As shown schematically in figure O a voltage-biased SPM 
tip is brought into close proximity with a graphene sample. Via capacitive coupling, the 
tip creates an image charge as a local change in the graphene charge density, indicated 
by the red spot, that adds to the average charge density n controlled by a back gate 
voltage Vg^ and to the existing disorder in the graphene layer, shown in blue, which 
consists of randomly placed scatterers created by charged impurities located above or 



below the graphene sheet [3J, l33]. As discussed below, the effect of the tip is to add one 
additional scatterer in the sample which can be moved about at will. 

A scanning electron micrograph of a graphene sample is shown in figure [3^. The 
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Figure 3. (a) Scanning electron micrograph of a graphene Hall bar, contacted by 
six Cr/Au leads, (b) Schematic diagrams of the graphene samples measured here (to 
scale). Sample C continues to the left and right for several microns. An ac current is 
applied between contacts 1 and 2 for all three samples, and the voltage is measured 
between contacts 3 and 4 for samples A and B, and between contacts 1 and 2 for 
sample C. 
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samples studied in these experiments are fabricated from single-layer graphene flakes 
deposited by mechanical exfoliation (the "sticky tape method" ) on a degenerately doped 
Si substrate capped with 280 nm of Si02. The samples were contacted by Cr/Au leads 
defined by electron beam lithography, and the graphene structures were formed by a 
mask defined by electron beam lithography followed by an oxygen plasma etch. A planar 
back gate voltage Vg is applied between the sample and the conducting Si substrate, 
allowing us to tune the charge density n in the graphene. 

The data presented here are from three graphene samples, shown schematically in 
figure [3b; all showed the same qualitative behavior. Samples A and B were measured in 
a four-probe geometry, with lead 2 grounded, and a root-mean-square (rms) current 
/ = 25 nA at 5 kHz applied between leads 1 and 2; the voltage for conductance 
measurements was measured between leads 3 and 4, with a lock-in amplifier. Sample C 
was tested using a two-probe geometry with the same current applied between leads 1 
and 2, with the voltage also measured across leads 1 and 2. For sample C, we estimate the 
contact resistance Rq by finding the value of Rq that fits the characteristic conductance 
expression G = 1/{R — Rq) vs. V^, which is linear on either side of the Dirac point 
yuirac] similar data is shown below in figure UK for the four-probe geometry. We find 
Rq = 7500 f^, which we subtract from the data from sample C. 



Each sample is mounted in a home-built scanning probe microscope j2J, 136|] , and 
cooled in He exchange gas that is in thermal contact with a liquid He bath at T = 4.2 K. 
The sample sits in the core of a superconducting magnet that provides a magnetic field 
up to 5 = 6 T perpendicular to the sample plane. We verify that the samples are 
single-layer graphene by observing quantum Hall conductance plateaux at the expected 
values of 4(z/ + l/2)e^//i, where u is an integer. For the scanning probe measurements, 
a conducting SPM tip with radius of curvature rup = 20 nm is held at a constant height 
htip = 10 nm above the graphene sample. The tip is grounded, so the tip charge is 
set by the contact potential between the graphene and the degenerately doped Si tip. 
Charged impurities on the surface of the graphene layer create image charges in the tip, 
which also contribute to the tip charge. 

Using transport measurements and electrostatic simulations, we characterized our 
graphene samples, and determined the spatial profile of the density perturbation created 
by the SPM tip. Classical electrostatic simulations (Maxwell, Ansoft LLC) were used 
to determine the charge density profile in the graphene created by the back gate 
and the SPM tip. The back gate voltage Vg capacitively creates a uniform charge 
density n = a{Vg — Voirac)) with a = 8 x 10^^ V~-^cm~^, where Voirac is the offset of 
the Dirac point from = caused by charged impurities near the graphene. The 
mobility ji for samples A, B, and C was jia = 7200 cm^/Vs, /hb = 5600 cm^/Vs, and 
fic = 4200 cm^/Vs, found from the measured G vs. Vg curves. The local perturbation 
to the graphene charge density caused by the tip has a Lorentzian-like shape with half- 
width at half maximum (HWHM) ^ 25 nm, and a maximum ^ 5 x 10^^ cm~^ for a 
potential diflFerence ^ 1 V between the tip and the sample. Previous experiments using 
scanning tunneling microscopy [35|, l37|] and a scanning single-electron-transistor charge 
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sensor [29| have observed the disorder in graphene samples on a Si02 substrate to consist 
of charge puddles with peak charge density ^ 4 x 10^^ cm~^ and lateral dimensions 
/ ^ 20 nm. The perturbation created by the tip in this paper has approximately the 
same amplitude as these naturally occurring charge density fluctuations, and about 
twice the width. Therefore, it is reasonable to think that the tip adds an additional, 
controllable scatterer to the pre-existing scatterers created by charged defects above and 
below the graphene. 




B (mT) B (mT) 

Figure 4. (a) Conductance G of graphene sample A measured in a four-probe 
geometry vs. back gate voltage Vg dit B = 0. Red line is a linear fit of G vs. Vg 
below 20 V. (b) Fluctuations AG vs. Vg about the linear fit {B = 0); the traces 
repeat when the SPM tip is at a fixed position far from the sample, (c) Same as 
(b), but with the tip at height hup = 10 nm above the sample for three different 
lateral positions spaced 100 nm apart, showing a clear dependence of the conductance 
fluctuations on tip position, (d) Magnetoconductance traces G vs. B with Vg = V 
(n = — 1.8 X 10^^ cm~^) and hup = 10 nm, at three tip positions spaced 5 nm apart; 
the traces repeat, (e) Same as (d), with four tip positions 100 nm apart; the traces 
now differ. 



4. Experimental Results 

Using an SPM tip, we probe the dependence of coherent magnet otransport on the 
position of a single scatterer. Figure H^i shows the conductance G of sample A vs. Vg 
at zero magnetic field; the data vary linearly with Vg on either side of the Dirac point 
VDirac = 22 V as indicated by the linear fit shown in red, as expected for single atomic 
layer graphene. The small fiuctuations AG away from the linear fit are reproducible, 
and are attributed to UCF [16]. The fiuctuations AG can be seen more clearly by 
subtracting the linear background from G, as is shown over a small range of Vg in 
figure Hb. The three traces of AG vs. Vg repeat the same scan with the SPM tip fixed 
far from the sample (tip height hup > 100 /xm), demonstrating good reproducibility. 
The rms amplitude 5G ^ e^//i is in agreement theory jsl, |9|. Figure Hb shows the 
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effect of moving the charged SPM tip to three different positions spaced 100 nm apart, 
at htip = 10 nm above the sample - the traces are now quite different. These data 
demonstrate that the motion of a single scatterer created by the SPM tip is sufficient 
to rearrange the conductance fluctuation pattern. 

Moving the tip also changes fluctuations in the magnetoconductance G. Figure [4]1 
shows three traces of G vs. B with l^ = 0(n = — 1.8x 10-^^ cm~^) and with the tip at 
three nearby positions, only 5 nm apart. Here the fluctuations repeat, with amplitude 
5G ^ /h. In contrast, when the tip is moved to several locations spaced 100 nm apart 
(figured^), the magnetoconductance is largely uncorrelated and looks like a different 
sample. Clearly, there is a critical length > 5 nm beyond which the tip must be moved 
in order to decorrelate the magnetoconductance fluctuations. 



4-1' Conductance images and correlations in a magnetic field 

We obtain conductance images G(5, Vg) by recording G as the tip is raster scanning 
over a 400 x 400 nm^ region in the center of the sample, at fixed magnetic field B and 
back gate voltage Vg. A conductance image G consists of 80 x 80 pixels that display 
the sample conductance G(r) = Gy^^ over a two-dimensional array of tip positions 
Yi^j = (x^, i/j) with pixel spacing Ax = = 5 nm, as shown schematically in figure [5b. 
The standard deviation SCup of the conductance G in an image is given by the standard 
deviation of Gr^j over all values of i and j. By acquiring a series of images G(5, Vg) as B 
or Vg is stepped to different values, we study the effect of the tip on coherent transport. 

A series of conductance images G{B^Vg) vs. B is shown in figure with Vg = 
(n = — 1.8 X lO"*^^ cm~^). The conductance images show fluctuations AG vs. r with 
standard deviation SCup ^ /h and lateral sizes ^ 10s of nm comparable to the Fermi 
wavelength and to the size of the scatterer created by the tip. The conductance images 
G(S, Vg) repeat over time intervals ^ 1 hr. as expected for coherent conductance 
fluctuations, demonstrating that their origin is not temporal noise. It can be seen by 
eye that the series of conductance images G(5, V^) shown in figure change with B 
in a continuous fashion. We quantify this change by calculating the cross-correlation 
Cab = J {Ga{^) — {Ga)){Gb{^) — {GB))dr between conductance images G^ = GA{^i,j) 
and G^ = GBi^ij)^ where angle brackets denote the average over the tip position r. 
We then define a normalized correlation Cab^ such that the autocorrelation of an image 
with itself is unity 

Pi _ Cab /r\ 

The grayscale map in figureEb displays correlations C(^Bi){B2) between images G(5i) 
and G(52) taken from a series of 200 images spanning B = ±145 mT at = (as 
in figure [5^). By definition, C(^b){b) = 1- The fainter peak along Bi = —B2 shows 
that there is some symmetry in the magnetoconductance fluctuations as the sign of 
B is reversed about B = 0. The symmetry is imperfect because of the nature of a 



coherent four-probe measurement [33|] , and because of slow drift of the images over the 
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Figure 5. (a) A series of conductance images G{B^Vg) vs. magnetic field B at fixed 
back gate voltage Vg = V (sample A), (b) Schematic of a conductance image G, 
with pixels spaced 5 nm apart displaying the conductance with the tip located 
at position Vij = (xi^i/j). (c) Correlation C(^b^){B2) between conductance images 
G(5i, Vg) and G(52, Vg) with = V. Lines show contours of a best fit from which 
/S.Br is obtained; the width of the black bar is 2ABr. 



course of the measurement time 10 hrs. We fit the map in figure Eb to a function 
given by the sum of two Gaussians centered at Bi = ±^2, shown as the contour plot 
in figure [5b. From this fit, we extract the magnetic correlation length ABc = 17 mT, 
seen in figure [Sb to be constant over the range of B shown here; ABc represents the 
magnetic field change needed to reduce the correlation between conductance images by 
one half. Using equation [U we find the diffusive coherence length L^^c = 500 nm, where 
L(fyc represents a value of L^^ obtained from ABc. A series of images G(5, Vg) vs. B 
provides a means to obtain ABc and L^j^c at a fixed back gate voltage V^, which we will 
return to below to study the behavior of ABc ^s. Vg. 

The amplitude SGup of conductance fiuctuations created by moving the local 
perturbation of the SPM tip differs from the amplitude SGvq created by changing the 
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the global back gate voltage Vg. This is understood by considering the ratio L^/L of the 
diffusive coherence length L^^ to the sample length L Q . A sample with width W ^ L^f^ 
and length L > L^j^ can be envisioned diS N = L/L^f) uncorrelated, coherent regions in 
series, each with conductance Gi = Go and exhibiting coherent conductance fluctuations 
with standard deviation 5Gi. The total conductance G of the sample is 

For a global change such as a change in V^, all regions undergo uncorrelated 
conductance fluctuations with standard deviation 6Gi = 6Go. The standard deviation 
SGvg in the total conductance is given by taking partial derivatives of equation [6] summed 
in quadrature: 



5G 



Va 



(7) 



Alternatively, the perturbation from the SPM tip acts locally, causing conductance 
fluctuations in only one coherent region yielding 5Gj = SGq and 6Gi = for i j . For 
a local perturbation, the standard deviation SGup in the total conductance is 

1/2 



^Gtip — 



E 



dG 
dG, 



SGi 



N-^5Go = {L^/LfSGo. 



(8) 



Combining equations [71 and [8] we obtain SGup/SGvg = {L^f^/Ly^'^. 

We obtain the standard deviations SGup and SGvg experimentally from conductance 
images G and traces of G vs. Vg for sample A. To reduce the uncertainty, we average 
SGtip over 20 conductance images with back gate voltage stepped from = — 1 to 1 V 
at 5 = yielding SGup = 0.32 e'^/h. We calculate SGvg = 0.50 e'^/h from the standard 
deviation of conductance traces G vs. Vg in the same range of = —1 to 1 V, measured 
with the tip away from the sample and with the linear background subtracted, as in 
flguredh. We flnd that the fluctuation SGup created by the tip is smaller than that from 
the back gate voltage by the ratio SGup/ 5Gvg = 0.64. This value is in good agreement 
with theory, which predicts SGup/ 5Gvg = [L^^jLY^^ = 0.6, using L^^ = 500 nm obtained 
above at K = 0. 



4^2. Weak localization 

Weak localization can be identifled by the conductance dip (or peak) that occurs at 
zero magnetic fleld. A conductance dip AGwl created by weak localization is shown 
in the magnetoconductance traces in flgure [61 As for other measurements of exfoliated 
graphene 10,14, 21], we observe a conductance dip from weak localization, not a peak. In 



our samples, where the sample size is comparable to L^^, the conductance change AGwl 
from weak localization is comparable to that from UCF. To observe the weak localization 
effect alone, we average over multiple magnetoconductance traces, by varying either the 
density n or the tip position r to average out UCF. 
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Figure 6. (a) Weak localization conductance dip in sample A measured with the 
SPM tip fixed far from the sample. Gray: Magnetoconductance AG vs. 5 at 11 
densities stepped from n = —1.6 x 10^^ cm~^ to n = —2.0 x 10^^ cm~^; curves are 
shifted to AG = at 5 = 0. Black curve: Average of the 11 gray curves. Red: 
Fit of black curve to a Lorentzian, yielding the width ABwl = 8.1 ± 0.5 mT. (b) 
Weak localization dip in sample A measured at tip height hup = 10 nm. Gray: 
Magnetoconductance G vs. 5 at a set of different tip positions separated by 100 nm, 
with fixed density n = —1.8 x 10^^ cm~^. Black curve: Average of 6400 such curves 
with tip positions spanning a 5 nm grid. Red: Fit of black curve to a Lorentzian, 
yielding width ABwl =9 ±2 mT. 

We clearly observe the weak localization dip at 5 = in figure [6K, by averaging 
magnetoconductance traces at different densities n with the tip away from the sample 
{htip > 100 fim). The gray lines in figure show a series of 11 magnetoconductance 
traces for densities n stepped by An = 4 x 10^^ cm~^ between n = — 1.6 x 10^^ cm~^ 
and n = —2.0 x 10^^ cm~^, a range AVg = 5 V. Because G changes significantly with 
n, these traces have been shifted to pass through the same point at 5 = 0. A weak 
localization dip AGwl is clearly seen in the average of all 11 traces, shown by the black 
curve, with magnetoconductance fluctuations averaged out. 
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We also observe the weak localization conductance dip AGwl in figure (6)3 by 
averaging magnetoconductance traces Gr^j vs. B over a grid of different tip positions 
Yij. The gray lines show magnetoconductance traces at fixed n = — 1.8 x 10^^ cm~^ and 
different tip positions, on a grid spaced by Ax = = 100 nm, with hup = 10 nm. This 
data is obtained from a series of 200 conductance images G(5, Vg) at different magnetic 
fields S, stepped from B = —145 mT to B = +145 mT. A magnetoconductance trace 
Gr,j{B) at tip position r^j is extracted from the dataset by taking the same pixel from 
each image G(5, Vg) at different values of B. The weak localization conductance dip 
AGwL is visible in the black curve in figure [6b, which is the average of all 6400 traces. 
Conductance fiuctuations are still clearly visible in the averaged magnetoconductance, 
because the tip-created scatterer only affects the conductance in a region of size ^ L^^, 
about one third of the sample. Magnetoconductance fiuctuations arising from regions 
of the sample at distances > L^^ from the tip are not decorrelated by changing the tip 
position. 

By measuring the width ABwl of the weak localization conductance dip AGwl 
shown in figures [6K and [6b, we obtain an estimate for the coherence length L^f^wL- 
To obtain ABwl^ we fit the averaged curves of G vs. B in figures [6^1 and [6b to a 
Lorentzian function, which approximates the predicted lineshapes given in equations [2] 
and [3] for weak localization in two dimensions 0, [l2|- The red lines in figures [6K and 
[6b show the results of the Lorentzian fit, with HWHM ABwl = 8.1 ± 0.5 mT and 
ABwL = 9 ± 2 mT in figures [6tt and b respectively, in good agreement with each other. 
Using equation [U we obtain a coherence length L^i^wl = 510 ± 20 nm for figure [6^ 
and L(i)WL = 480 ± 50 nm for figure [6b. These values obtained from weak localization 
measurements match the estimate L^^c = 500 nm obtained above from cross-correlations 
between magnetoconductance images shown in figure [5l 

Measurements of ABwl from the weak localization conductance dip AGwl and 
measurements of ABc from cross-correlations between conductance images allow us to 
study how the phase coherence length L^^ depends on the charge density n. Previous 
bulk weak localization measurements found the coherence length L^wl in graphene 



varies significantly with n, with a minimum at the Dirac point |j, |2lt]. Figure [7K shows 
our data for ABwl and ABc/2 vs. n, obtained from conductance images G(5,V^) 
for sample C. The agreement clearly shows the factor of 2 difference, predicted by the 
interference of two time-reversed diffusive loops for weak localization vs. the single loop 
for interfering forward scattering paths, discussed above. Each data point is derived 
from a series of 200 conductance images G(S, Vg) at a fixed density n. The circles show 
ABwl-, from the width of the weak localization conductance dip in an average over tip 
positions, as in figure [6b. The squares show ABc/2 obtained from the cross-correlation 
C(Bi)(B2) (equation [5j) between conductance images G(Si,V^) and 6(^2, V^). We find 
ABc from the correlation function C(^Bi)(B2) as in figure [5t. The values of ABwl and 
ABc in figure [7K both increase as the hole density decreases and the system approaches 
the Dirac point. 

Figure [7b shows ABwl found for sample B with the tip fixed far from the sample. 
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Figure 7. (a) Widths of the weak locahzation conductance dip ABwl and the 
correlation magnetic field ABc for conductance fiuctuations vs. charge density n in 
sample C. Circles: ABwl vs. n obtained from magnetoconductance traces averaged 
over different tip positions, as in figure [fib. Squares: ABc/2 obtained from cross- 
correlations between conductance images, as in figure [5b. (b) ABwl in sample B, 
obtained with the tip fixed far from the sample, by averaging magnetoconductance 
traces over a range (4 x 10^^ cm~^) of n, showing qualitative agreement with the data 
in (a), (c) Standard deviation SGup of the conductance Gr^j over all tip positions r^j 
in a conductance image G vs. n for sample A. 



Each data point is obtained by averaging 20 magnetoconductance traces G vs. B over a 
range of densities n = ±2 x IQ-^-^ cm~^ (a 5V change in Vg). The data in figure [7b show 
that ABwl increases as the electron density is decreased, and is largest at the Dirac 
point (n = 0), complementing the behavior for holes. 

The coherence length for samples B and C is estimated from ABwl and ABc 
in figures [TK and [Zb, using equations [H and [H We find that the minimum coherence 
lengths LcfyWL = L^f^c = 300 nm occur at the Dirac point for both samples. The maximum 
coherence length occurs at the largest hole or electron density: L^i^wl = 600 nm for 
n = 4 X 10^^ cm~^ in sample B, and L^f^wL L^f^c 500 nm at n = —1.7 x 10^^ cm~^ in 
sample C. The change of L^p with n is in good agreement with previous results (4l. l2l[| . and 
is likely caused by the change in electron-electron interaction strength with n jj, ll7|, [21 



22 



as well as the effects of electron-hole puddles near the Dirac point 

The standard deviation SGup of conductance fluctuations over a single image for 
a series of conductance images G{n) at different densities n at 5 = is shown in 
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figure [Tt. We find that SGup increases away from the Dirac point with the electron 
or hole density; these results are similar to bulk transport measurements on single and 
multilayer graphene 19, 22, 38|. The trend in SGup vs. n is consistent with the behavior 
of ABwL and ABc in figures [TK andlZb- As discussed above, when the sample length 
L > L(f)^ then SGup oc {L^f^/Ly. Because ABwl oc L^^^ and ABc oc L^^^ we expect 
6Gtip oc 1/ ABwL and SGup oc 1/ ABc] the data in figure [71 follow this trend. From n = 
to n = ±2 X 10-^^ cm~^, SGup changes by a factor ^ 2, which agrees with the observed 
change in ABwl and ABc by a factor ^1/2 over the same range of n in figures [7K and 

03. 



5. Conclusions 



The study of coherent electron transport and weak localization in graphene has recently 
received much attention through bulk transport measurements. Here, we add a new 
technique that allows one to probe transport by creating a movable scatterer with an 
SPM tip. This allows us to investigate weak localization in a new way by creating 
a controllable, local change to the disorder configuration. Our measurements provide 
a direct way to observe the spatial nature of coherent transport in graphene. This 
technique provides an image of the magnetoconductance vs. tip position that represents 
a spatial "fingerprint" of the interfering paths at a particular Fermi energy and 
magnetic field B. From correlations between images at different B we measure the 
coherence length L^^c for electrons diffusing through the sample, and find good agreement 
with the coherence length L^j^wl obtained from weak localization measurements. 

Understanding and controlling disorder is one of the main challenges in realizing 
many of the proposals for probing new physics and discovering applications for graphene. 
In this work, we have described a tool that maps the effect of a nanoscale change to 
the disorder in a graphene sample produced by an SPM tip. This represents a step 
towards the goal of gaining mastery over disorder in graphene, by developing methods 
to controUably shape it and bend it to our will. 
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